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The present work studies the Ghatak-Sherrington (GS) model in the presence of a magnetic 
random field (RF). Previous results obtained from GS model without RF suggest that disorder 
and frustration are the key ingredients to produce spontaneous inverse freezing (IF). However, in 
this model, the effects of disorder and frustration always appear combined. In that sense, the 
introduction of RF allows us to study the IF under the effects of a disorder which is not a source 
of frustration. The problem is solved within the one step replica symmetry approximation. The 
results show that the first order transition between the spin glass and the paramagnetic phases, 
which is related to the IF for a certain range of crystal field D, is gradually suppressed when the 
RF is increased. 
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I. INTRODUCTION 

In recent years, there is considerable interest in an un- 
usual class of phase transitions known as Inverse Tran- 
sitions (IT). In these counterintuitive transitions, the 
ordered phase has more entropy than the disordered 
one. It is important to remark that such phase transi- 
tions have been observed experimentally in physical sys- 
tems as, for example, magnetic thin films 1 and high-Tc 
superconductors^. From the theoretical point of view, 
the understanding of the puzzling mechanisms behind 
such transitions is a challenging task. In this sense, some 
magnetic models can be quite useful (see, for instance, 
Refs. 1 3J-TlXr> - Indeed, in some of these mentioned refer- 
ences, it is suggested that disorder and frustration are 
key ingredients to produce spontaneously IT and, in par- 
ticular, Inverse Freezing (IF). In this kind of IT, the spin 
glass (SG) phase is the ordered one. Nevertheless, in 
most of these spin models, disorder appears only as a 
random bond with no other forms of disorder being con- 
sidered. 

A new form of disorder that could be added to these 
models is a random field. Albeit the addition of a ran- 
dom field in these spin models is a new simple term, it can 
introduce a complex interplay of two highly non-trivial 
manifestations of disordered magnetism, i. e., spin glass 
and random field s 12 ' 13 . As an example of this complex 
interplay, recent studies of the Shcrrington-Kirkpatrick 
(SK) model^ with a random field indicate that the SG 
phase is deeply affected by the presence of such field. 
In particular, not only the nature of the SG boundary 
phase can be affected^— but also its location can be 
change d 18 ' 19 . On the other hand, it should be noticed 
that to increase the disorder in the problem by adding 



a random field docs not mean necessarily to increase 
frustration. In fact, as shown by the Random Field 
Ising model (RFIM) in the mean field approximation, 
the presence of a random field is an additional form of 
disorder, nonetheless, it is not an additional source of 
frustratio n 20 ' 21 . 

An example of a simple disordered magnetic model 
presenting spontaneously IF in the mean field approx- 
imation is the Ghatak-Sherrington (GS) model 2 ^. The 
precise meaning of spontaneous existence of IF is that in 
this model it is not used any artificial procedure to in- 
crease the entropy of the SG phased. From the previously 
mentioned results for the SK model and the RFIM, one 
could ask what would happen with the spontaneous IF 
found in the GS model when a random field is applied? 
Is the reentrance in the first order boundary between the 
SG and paramagnetic phases affected by the presence of 
the random field? The answer to these questions can not 
only be useful to clarify the role of disorder as an ingredi- 
ent to produce (or destroy) an IF, but also to distinguish 
the effects coming from disorder, from those ones com- 
ing from frustration. The reason is that, at mean field 
level, the effects of a local disorder coming from random 
fields can be distinct from those ones coming from ran- 
dom bonds. For this last case, the effects of disorder and 
frustration can appear combined, as in the GS model. 

In the present work we study the GS model in the 
presence of a longitudinal random field. In this model, 
the spins can assume values S = ±1,0, D is a crystal 
lattice field and the spin-spin coupling J^ and the ran- 
dom field hi follow Gaussian distributions. In the case 
of a random field following a Gaussian distribution, it is 
necessary to use the replica symmetry breaking in order 
to obtain the SG solutionis. The free energy is obtained 



and consequently, the SG order parameters. The disor- 
der is treated using the replica method in which the SG 
solutions are searched within the one-step replica symme- 
try breaking (1S-RSB) approximation^. The location of 
the first order SG/PM boundary phase, in the present 
work, is done by using the same criteria proposed in Ref. 
[24 . In this case, the chosen SG solution in the transition 
is that one which meets continuously the SG solution at 
small D. This solution corresponds to the lower SG free 
energy. It is important to point out that the GS model, 
besides disorder and frustration, has the ingredients, as 
proposed by Schupper and Shnerb^, to present sponta- 
neous IT, i.e., an interacting term favoring a magnetic 
state and a mechanism favoring non-magnetic states in 
the sites. 

It should be remarked that in the GS model there 
is a new order parameter p, which is related to the 
diagonal replica components (the spin self-interaction 
p = ((S a S a )) j : h.) ■ In fact, this order parameter is re- 
lated to the average number no = 1 — p of non-magnetic 
states in the sites which also reflects the capability of 
the sites to interact or not2£. For example, the limits 
no = and no — 1 can be related to the fully interacting 
(p = 1) and the fully noninteracting (p = 0) regimes, 
respectively. The behavior of no can be important to 
clarify the mechanism behind the existence of the IF in 
the present approach. 

The paper is structured as follows: in section II, the 
free energy within 1S-RSB scheme was found. In section 
III, a detailed discussion of phase diagrams is presented. 
The last section is reserved to the conclusions. 



II. GENERAL FORMULATION 



The model is given by the Hamiltonian below 
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where the spin can assume values S = ±1,0, the summa- 
tion (i,j) is over any two interacting sites and D is the 
crystal field. The spin-spin coupling Jy and the mag- 
netic fields hi arc random variables following independent 
Gaussian distributions 
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Herein the procedure introduced in the Refs. [TafTg is 
closely followed. Therefore, the free energy per spin can 
be obtained using the replica method as 

-f3f= lim 1im-±-(«Z n »j, h -l), (4) 

AT->oo n-¥0 Nn 



where the average << /({.%} , {hi}) »j,h is given by 

« f({j tj } , {hi}) »j, h = f n [dj l3 p{j l3 )] 
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with Jij and hi following the Gaussian distributions de- 
fined in Eq. [I] 

From Eq. [1J after averaging Z n over J^ and hi, one 
gets by using the Hubbard-Stratonovich transformation 
the following relation: 



H&f. 



(/3J) 2 



(G) 



^ ? L+lnTre L , 

a ) 

where L is the effective Hamiltonian 
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and the diagonal q aa and non diagonal q a p elements of 
the replica order parameter matrix are given by 



q aa — < S a >, q a (s —< S a Sp >, 



(8) 



respectively. The symbol < ... > means the thermody- 
namical average over the effective Hamiltonian L given 
in Eq. 

The Parisi's 1S-RSB^ is used in Eq. H Therefore, the 
replica matrix q is parametrized as shown below 2 ^ 



q a p = 9i if I(a/x) = I(/3/x) 



q a p = qo if I(a/x) ^ I(/3/x) 
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where I(x) gives the smallest integer which is greater 
than or equal to x. While 



q aa = p. 



(11) 



The parametrization introduced in Eqs. l9lllll is used in 
the sum over the replica labels in Eqs [6][7] 

The final step in the 1S-RSB is a linearisation in Eq. 
[7] which introduces new auxiliary fields z and v. Finally, 
the free energy is given by 
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where 

K(z, v) = 1 + 2e 7 cosh h(z, v) 
with f Dw = f J™ er v> 2 /'* ( w = z v ) 
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and h(z,v) is 
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The equations for the order parameters go, Qi, P and 
the block size parameter x are obtained from Eq. [12] 
by using the saddle point conditions. Other thermody- 
namic quantities can also be derived from Eq. [12l as for 
instance, the entropy s = —df/dT. 



III. RESULTS 

In this section, the numerical results for the order pa- 
rameters are displayed with different set of values of the 




FIG. 1: Order parameters as a function of temperature for 
D/J — (upper panel), D/J = 0.5 (lower panel) and several 
values of A. 

width A (see Eq. \S§ and the crystal field D that are 



given in units of J. These results allow us to build sev- 
eral phase diagrams in which the SG phase is character- 
ized by the RSB order parameter S = q\ — go > 0. For 
A > 0, the paramagnetic (PM) phase is given by S = 
but gi = qo =g^0. The behavior of entropy and ng^ 
are also presented. Figure Q] exhibits the behavior of or- 
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FIG. 2: Phase diagrams of temperature versus A/J for 
D/J = 0.0 (full line) and D/J = 0.5 (dashed line). 

der parameters for different configurations of A and D. 
In particular, the values chosen for D corresponds to a 
range in which a PM/SG second order phase transition 
appears. For A = 0, it is obtained, as well known, that 
the temperature where 5 > (the freezing temperature 
Tf) is decreased when D increases^. The results also 
indicate that p depends on the temperature, with every 
site becoming magnetic asT->0 for D = (j> —> 1). 
Most important, for D > 0, p < 1 as T — >• 0, which 
indicates that D affects directly the average number of 
non-magnetic sites no. When D = 0, the results are sim- 
ilar to those ones found recently for the SK model with a 
random field (see RefM-). For A = 0, the PM/SG phase 
transition is trivially induced. However, for A ^ 0, the 
random field plays a role similar to a uniform magnetic 
field suppressing the SG solution in the replica symmetry 
(RS) scheme. In that case, the PM/SG phase transition 
can only be found in the 1S-RSB scheme (6 > 0). One im- 
portant consequence is that Tf decreases as A enhances. 
However, the spin self-interaction p remains dependent 
on temperature but much less affected. Particularly, at 
very low temperatures. Finally, when the effects of finite 
D and A are combined, some of the previous effects are 
preserved but others are more pronounced. For instance, 
the random field still acts as a uniform magnetic field 
suppressing the SG solution. Despite that, Tf decreases 
faster and p is more affected with the increase of A for 
intermediated temperatures than the case D = 0. 

The behavior of Tf discussed above is better analyzed 
in the phase diagrams of Fig. [2] For instance, Tf is 
gradually decreased by A. However, the SG phase is al- 
ways found at lower temperatures (see full line in Fig. ^ . 
For D/J = 0.5, Tf appears at even lower temperatures. 
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FIG. 3: Phase diagrams of temperature versus D/J for several values of A. Panels (b), (c) and (d) show details of the first 
order transitions. The dashed-dotted lines at low T are extrapolations to the T = 0. 
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FIG. 4: Entropy as a function of temperature (region where 
occurs the reentrance) for three values of A and D. The first 
order temperature is equal to Tif = 0.1J for all curves. 



However, the results suggest that the SG phase can still 
be found even for very large values of A. 

Fig. [3] shows the results for larger values of D. In 
that case, the continuous PM/SG phase transition is de- 
creased until a tricritical point as D enhances. After the 
tricritical point, the phase transition becomes first order 
where there are more than one PM and SG solutions. 
The first order boundary is then located by comparing 
the free energy of the stable PM solution to the 1S-RSB 
free energy of the SG solution that presents lower free 
cncrgy24 Particularly, the SG/PM first order transition 
can be reentrant. This reentrance represents an IF, in 
which the usual entropic relation between the ordered 
and disordered phases is changed (see Fig. 2]). This re- 
sult for A/ J = 0.0 was already discussed in Ref. [j. On 
the other hand, the presence of the random field affects 
significantly the first order phase transition line. The 
tricritical point is displaced to lower temperatures and 



higher D as A increases. For A large enough, the tri- 
critical point is not found in our range of temperatures 
as well as the SG/PM reentrance. However, for smaller 
values of A, the reentrance is still observed for a certain 
range of D (see Figs. Efb)-(d)). Nevertheless, the reen- 
trance is gradually diminished by A. It suggests that 
although the disorder introduced by the random field is 
relevant to the location of critical lines, it is not able to 
enhance the IF region. On the contrary, the reentrant 
transition tends to be suppressed by the presence of a 
random field. 

Fig. 2] presents the entropy s as a function of T/J for 
three values of A/ J: 0.0, 0.25 and 0.40. The values of 
D/J are chosen in an order that they allow us to follow 
the entropy difference between the SG and PM phases at 
the same first order boundary temperature T\f = 0.1 J. 
For A = 0, the entropy of the PM phase is found below 
the SG one at the first order transition. As already pre- 
sented in Ref. y, this result confirms the existence of IF 
in the GS model. When there is a random field, the en- 
tropy at low T/J increases especially the PM phase. As 
a consequence, the entropy difference of the SG and PM 
phases at the transition is decreased. This result con- 
firms that the reentrance seems to be weakened by the 
random field. 

Fig. exhibits hq which can be useful to understand 
the behavior of phase diagrams T/J versus D/J shown 
in Fig. H For the isotherm T/J = 0.1 with A = and 
0.25 shown in Fig. Ufa), no increases continuously with 
D until the PM/SG first order phase transition where 
it tends discontinuously to the higher uq limit (see full 
line in Fig. Ufa)). The behavior of no in the SG is not 
affected by the random field. However, for A = 0.25, 
the first order phase transition is displaced for a larger 
value of D as shown in Figs. GHb)-(c). For A = 0.5, the 
tricritical point is located below the isotherm T/J = 0.1. 
Therefore, the PM/SG phase transition becomes a sec- 
ond order one, with no increasing continuously. However, 
no in the PM phase is decreased as compared to the pre- 
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FIG. 5: Panel (a) shows the average number of sites with nonmagnetic states no versus D/J for T/J = 0.1 and several values 
of A/ J. Panel (b) shows no versus T/J for three values of D/J and A/ J. The values of D/J in panel (b) are chosen in order 
to present the first order temperature equal to Tij = 0.1J. Tf represents the freezing temperature. 



vious cases. This behavior could explain the presence of 
the SG order at a large range of D. To put in another 
way, D favors the nonmagnetic states (no — ^ 1), but the 
random field enforces the magnetic ones. Therefore, the 
random field can diminish the effects caused by D. 

Fig. [5][b) shows that no also depends on the tempera- 
ture. For A/ J = 0, no in the PM phase increases as the 
temperature diminishes until the Tf , where the SG order 
appears. However, the PM phase can be found again at 
lower temperatures. Most important, the behavior of n 
in the PM phase at low and high T are different. At low 
T, the no of the PM phase presents high level of nonin- 
teracting states as compared to its counterpart at high T. 
As a consequence, the reentrance can appear. In other 
words, the nontrivial frustration and the presence of non 
interacting states can introduce an IF transition. On 
other hand, the disorder caused by the random field al- 
ters the no distribution in such way that inverse freezing 
is weakened. It tends to decrease the difference between 
the no of the PM and SG phases at the transition. In 
addition, the no of PM phase at low temperature is de- 
creased by A, which results in the increase of the entropy 
of the PM phase at low T/J, as shown at Fig. 2J 



IV. CONCLUSIONS 

In the present work, the Ghatak-Sherrington (GS) spin 
glass model has been studied with the addition of a ran- 
dom longitudinal magnetic field hi and a crystal field D. 
This study was carried out within the 1S-RSB scheme 
using a mean field approximation. It is known that spon- 
taneous Inverse Freezing (IF) can appear in this model 
depending on the value given for D^. The main purpose 
of the present work is to investigate the effects of the 
random field following a Gaussian distribution in the IF. 
Since, it does not introduce any additional frustration in 



the problem (at least, at mean field leve£) its presence 
can be helpful to elucidate what is the role of disorder as 
an ingredient to produce (or suppress) the IF as distinct 
of frustration. 

Our results show that the random field depresses the 
freezing temperature Tf. Moreover and most impor- 
tantly for the purposes of the present work, the random 
field decreases the reentrance in the first order phase 
transition line PM/SG. The mechanism responsible for 
such effect is concerned with the average number of sites 
with nonmagnetic states in the SG and PM phases. It 
is important to remark that in the absence of a random 
field, the PM phase at low temperature is characterized 
by higher number of sites with nonmagnetic states as 
compared with the PM phase at higher temperature. By 
contrast, the random field tends to diminish the num- 
ber of sites with nonmagnetic states of the PM phase at 
lower temperature. In that way, the SG phase is again 
energetically favored at lower temperature. 

Therefore, one can conclude that the disorder given 
by the presence of a random field does not enforce the 
IF. On the contrary, the presence of the random field 
tends to decrease the difference of the entropy of the SG 
and PM phases at the first order phase transition. As a 
consequence, the random field not only weakens the reen- 
trance, but also tends to suppress the PM/SG first order 
transition corresponding to the IF by lowering the tricrit- 
ical point. Although, our result is limited to a specific 
magnetic model, it could be an indication that frustra- 
tion more than the disorder, would be a key ingredient to 
produce such counterintuitive type of phase transition. 
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